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ABSTRACT 

Large amplitude oscillation of cosmic field that may occur right after infiation 
and in the decay process of weakly interacting fields gives rise to violent particle 
production via the parametric resonance. In the large amphtude limit the problem of 
back reaction against the field oscillation is solved and the energy spectrum of created 
particles is determined in a semi-classical approximation. For large enough coupling 
or large enough amplitude the resulting energy spectrum is broadly distributed, 
implying larger production of high energy particles than what a simple estimate of 
the reheating temperature due to the Born formula would suggest. 
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It is a common belief that coherent field oscillation gives rise to particle produc- 
tion coupled to the field, and ultimately the field is dissipated away if no further 
supply of field energy is provided. The phenomenon may occur in a number of co- 
somological situations, notably in the reheating stage of infiation and in the decay 
process of weakly interacting spinless fields (WISF) such as the Polonyi field in some 
supergravity models. Recently it has been recognized [|||] @] that presence 
of the parametric resonance in the corresponding classical system much enhances the 
field decay, in typical cases of the strong coupling or the large amplitude oscillation 
the decay proceeding much more rapidly than what is expected by the Born term. 
The effect is important in many ways, and may change our understanding of the 
exit problem of inflation and the WISF problem. For instance, the inflaton decay 
may rapidly end with non-thermal particle creation and the energy redistribution 
towards the thermal one may take place by much slower processes of interaction 
among produced particles. The new decay mechanism may thus drastically change 
previous estimate of the reheating temperature. 

In this paper we investigate the back reaction of particle production against the 
fleld oscillation and calculate the energy spectrum of created particles before effect 
of particle interaction becomes important. We do this by fully considering effect of 
cosmological expansion. Throughout this work the oscillating fleld is treated as a 
classical variable and the back reaction is incorporated by a semi-classical approxi- 
mation. The formalism of Ref . based on the Schrodinger picture of the quantum 
system coupled to the fleld oscillation enables one to solve the problem of both the 
back reaction and the particle spectrum in a compact way. Moreover, we are able 
to clarify in what sense one may effectively regard a pure quantum state as a mixed 
state, even if no interaction among produced particles is explicitly taken into account. 

Inflaton or WISF fleld is generically denoted in this paper by ^(t). It is assumed 
that after some cosmic time the fleld undergoes a phase of the damped oscillation 
around a stable minimum taken to be ^ = here, like ^(t) oc a~^^'^{t) sin {m^{t — to)) , 
with the mass of ^ fleld, and a{t) the cosmic scale factor. The ^ fleld has coupling 
to various matter and radiation flelds, and we consider a typical Yukawa coupling 
of massless quantum fleld ip of the form ^ /i^y^^ . The generic quantum fleld ip may 
have their own interaction, and the state of the universe may be in a thermal equi- 
librium if the interaction rate is larger than the Hubble expansion rate. What we 
consider in the rest of discussion, however, typically occurs in time scales shorter 
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than 0(10 — 100) x In this circumstance one may envisage that the quantum 

field (f obeys a free field equation with the periodic mass term /i^ (t) in the expanding 
universe. This should be a reasonable starting point of discussion if one considers a 
time duration much larger than the oscillation time of order 1 / and shorter than 
the interaction time, unless the interaction rate is much enhanced by enormously 
large number density of created particles. 

Behavior of the if system is unambiguously described by the density matrix, which 
is then decomposed as a direct product of independent spatial Fourier component 
p{k) due to the assumed spatial translational invariance and neglect of interaction. 
The bulk of the universe may be in a thermal equilibrium, but our main interest 
lies in the state of newly created (f particles by the ^ oscillation. Hence we may 
concentrate on this part of the density matrix, which is determined by the wave 
function of a quantum state. 

The wave function of our problem has a Gaussian form of exp[z^ q'i]/^/u . The 
function u{t) here obeys the classical equation of harmonic oscillator with u;^ + 
^^{t) ,LU = \k\ as a time dependent frequency squared. The initial condition to this 
solution is specified by a choice of the corresponding initial quantum state, which 
we assume to be the ground state with the coupling /i^ term switched off at t = to, 
hence M(to) = (^^/^) ""^^^ > ''^(to)/M(to) = • Take as a convenient basis the Fock 
space that diagonalizes ip particle number when the coupling term fiC, is switched 
off. Non-vanishing density matrix element in each Fourier component is then given 
1 by 
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The time dependent real parameters, K(t) and 6{t), are related to u(t). In the 
small K limit most important in application (3? and 3 indicating the real and the 
imaginary parts), 

^=^7TW '""^=— ' withD=(vr---H)/H . 
For a while we ignore effect of cosmo logical redshift, focussing on the short time 
behavior compared to the Hubble time. With the dimensionless parameters defined 
hj h = ^ and 6* = ^ , the classical oscillator equation for u becomes of Mathieu 



type 



— + (h-2ecos(2z))u = 0, (2) 
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with z = . The amphtude ^, hence 6 is taken constant in time for this part 

of discussion. As is well known, the classical solution in the parametric resonance 
region is asymptotically written as u{t) = e'^"^?*/^ P(t) , with A real and positive, and 
P{t) a periodic function. With this form one may use for parameters of the density 
matrix 

Attuj 4 

The produced number of particles given by J2i 2lp2i,2i is ^ l/K oc |mP , and 
exponentially grows in the parametric resonance region after performing time average 
over the quasi-period. On the other hand, the phase factor 6 that appears in off- 
diagonal elements of the density matrix is almost periodic with period of Air/m^ . 
For a sizable \6\ ^ 0[1] the parametric resonance readily occurs in a large parameter 
region of {6 , h) space except in narrow bands of the stability region. We first wish 
to improve the analytic formula for large 6 and derive a form of the density matrix 
useful for application. 

One can envisage the following analogue quantum mechanical model in order to 
estimate the growth factor A. The classical equation Eq.|^ is formally equivalent 
to the potential problem in one dimensional quantum mechanics, with the potential 
given by 6 cos(22;) and with the energy h/2. For h <C 16*1 there is an infinite succession 
of large potential barriers. The difference from the usual eigenvalue problem is that 
one allows unbounded solution, since one assumes here a boundary condition at some 
z = zo instead of the boundedness condition at ±oo. Rather than considering an 
infinite series of barriers at once, one may relate input amplitude to output amplitude 
after passing just one barrier and repeat this process n times, with n = [m^t/(47r)] 
a large integer. 

We thus introduce a transfer matrix T that connects the input amplitude to the 
left of one barrier (at —2:0) to the output amplitude to the right of this barrier (at 
zq) as follows; 

y Ur{±Zo) Uji{±Zo) J 

taking left- and right-movers as the basis of outgoing wave, zq is assumed positive 
and large. S(±zo) is a 2 x 2 matrix that transforms from two independent modes 
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i'fiz), ^pf{z), with the properties of ^p^{±Zo) = l,ip'^{±ZQ) = 0, ^ljf{±Zo) = 

, 'ijj'j^{±Zo) = 1 , to the left and right moving modes ul{z) , uji{z) ai z = ±zq. The 
transfer matrix can be diagonahzed by a similarity transformation, T = S^^T^S , 
giving diagonal elements t± with > 1 and = < 1 , both real for an 
unbounded solution. One can then clearly identify the growth factor by A = . 

Our approximation for large 6 consists of replacing the periodic barrier by a piece- 
wise inverted harmonic oscillators: thus for an interval of \z\ < | , cos(2z) — > 

1 — 2z'^ . The approximation gets better as \9\ increases. Let e = ^ — 1 = 

2 

^ — 1 . The left- and right-moving mode functions with the asymptotic behavior 
of (2v^2)(-i^^^^)/2 exp(^iv^ z^) are written in terms of the standard parabolic 
cylinder function Di,{z) as 



ul 



Despite of the finite range Zq = 7i/2 the argument of these functions contain a large 
factor 9^^^, hence one may use the asymptotic form of these functions. By using 
analytic properties of the standard parabolic cylinder function 0, one derives 

S(-.o) =( ^ " U ] , (6) 

\a* f3* J \ ur{zo) -u'jiizo) J 

a = ^ ^ ^, (3 = ^e-^'^/\ (7) 




The transfer matrix is further simplified, by introducing a constant x by u'ji{zq) 
XUr{zo) and using ul = u\. We obtain 

with e^**^ = ur[zq) /u\[zq) . Eigenvalues of the matrix are computed to give \t±\ = 
I |trT| ± ^ (tr — 1 . The eigenvector corresponding to the growing mode is 

of great interest. In the basis of left-right movers the eigenvector Ti{zq)v^ has a 
form of (y ,V*)'^ . The growing mode is thus an equal mixture of the left and 
the right movers like —V*ul + Vur, which becomes a single sinusoidal function 
asymptotically. Profound implication of this fact will be discussed shortly. 

The similarity transformation S that diagonalizes T may be written in terms of 
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the eigenvectors v±. A large number of repeated barrier crossover leads to 
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(10) 

assuming a purely imaginary x relevant later. The initial state corresponding to the 
quantum ground state is given by Vq oc (1 , iu)'^ , u = ^ , in this basis. The norm 
of the state vector is best determined in the basis of the left and the right movers 
that have the unit flux. Comparing the two norms, one finds that after many barrier 
crossing the initial state vector vq is amplified by 



1 jx^ + co')ix^ + \xn 

for X purely imaginary. Besides the usual factor jt+l" = e"^°'*+', it contains a pref- 
actor. 

One can simplify relevant formulas in the large 6 limit. First, the constant x is 
of order y/O and is purely imaginary: x = 

iVe {2zo + ^) . The factor X in Eq.[l| 

is large and of order VO: 

V n 29' J l + tan(^/'/2) 

From this the extra prefactor in the amplification formula Eq.[Tl| turns out of order 
unity to leading order of 6 and has a complicated energy dependence. In the rest 
of numerical analysis we neglect this energy dependent prefactor for simplicity. For 
e < of our interest the growth factor A ~ In |23?(a;e^*'^)|/7r is computed as 

1/7T/ ln21,, , 

A = ^ (1 - ^) + — + - In I cos V^l , (12) 

^ = ^ + 2-'v/^-^M^)- (13) 

In the spirit of large Q approximation we ignore, in what follows, the constant term 
In 2/7r, and the complicated logarithmic term, which should be a good approximation 
except in narrow stability bands of | cos^/)! < e"'^^''^'/^ . The instability region with 
e > is very much limited and we ignore contribution from this region. Thus the 
main result in this part of calculation is summarized as A = ^ (1 — • 
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The result of the growing mode T,{zq)v^ has an important consequence on behav- 
ior of off-diagonal density matrix elements. From the formula Eq.^ and the fact that 
the asymptotic form of the growing mode is given by a single sinusoidal function, it is 
found that the periodic part of the classical solution P{t) has zeros, hence the factor 
tan (5 oscillates between ±00. The oscillation period is small, given by 0[l/m^]. The 
time average over a larger time interval then gives vanishing off-diagonal elements 
since p2« ,2m oc e"^**^*^'"™) . Furthermore, neglected interaction among produced parti- 
cles is expected to readily break the delicate phase coherence present in the quantum 
state of our consideration. Thus even if the quantum (f system is in a pure state 
in the strict sense, it is reasonable to suppose that the system is effectively in a 
mixed state, described solely by the diagonal part of the density matrix. From this 
diagonal density matrix one may compute necessary quantities such as the average 
energy density, the energy distribution function, and more importantly an effective 
entropy which will be discussed later. For the rest of discussion we thus use the 
reduced density matrix with off-diagonal elements projected out. 

We now consider the cosmological situation in which the amplitude ^ may vary 
in time. In order to discuss the back reaction against the ^ field oscillation, it 
is necessary to consider both the cosmological redshift and the dissipation caused 
by particle production. For this purpose we extend the previous analysis, first by 
taking an infinitesimally small time period during which cosmological expansion may 
be ignored and previous formulas under periodic perturbation can be applied, and 
then by incorporating adiabatic change of the amplitude ^(t) within time integral. 
Thus we replace the previous factor like A(^)t by //^ dt' X {^{t')) . 

The ^ field equation with the back reaction included is derived by considering en- 
ergy balance between the ^ field and created radiation. Combined with the Einstein 
equation, it leads to a set of equations among the energy densities, p^(t) = \ m'^^'^it), 
Prit), and the scale factor a(t); 



where (p) = Tr (loN^pd) , with the trace including the phase space integration. 




(14) 



(15) 



(16) 



Because A = y ^ (1 — ^)/m^ , the source term ^ (p) for a single mode has the 
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form of 

- 31' ^''^'^ X* - 31 Si' " • 

with u!'^(t)a'^{t) independent of time. Somewhat by accident, the phase space integral 
involved here can be performed explicitly in the form of J^^ (icjcj^"+^ e^'^'^^ with 
n — 1,2. This yields a closed set of differential equations instead of complicated 
integro-differential equations. In terms of dimensionless quantities divided by the 
initial ^ energy density = | ^|Co > ^ = Pc/ p\-i = PrI p\-i = "^^^ ? these are 
(' indicating ^) 

Y'W n' n' 
U'^c^oY'/'U, W'^c^oY'/' + -^ + 2-W, - ^ dWTZ , (17) 

^1 a a 

Y' + 3dVY + ZY ^ -BY^'"^ U[W + {W + 2) ]/W^ , (18) 
Z' + Ad^/YTZZ = B Y^l^ [/ [ - 2 + e-^ (ly + 2) \IW^ , (19) 

where = ^ = ^o/2 . d ^ ^ B ^ {^)\ with c a numerical factor 

of order unity: c ~ l/\/2 . By rescaling the time variable, one sees that this system 
of equations has two intrinsic parameters, c^go/dP and B/d, which we call initial 
amplitude and damping factor in what follows. 

We first discuss asymptotic behavior of solution that can be worked out by ana- 
lytic means. A consistency check shows that the ^ and the radiation energy density 
asymptotically behaves like 

Pi ^ p^(c^Tln(T/Te))-^ Pr ^ P?/(2dT)^ a' / a ^ 1/(2t), 

where t^. is a constant determined from early time behavior. The created radiation 
density pr asymptotically follows the standard time variation of radiation dominated 
universe, but variation of the ^ energy density is balanced against both the redshift 
and the particle creation terms. Unless some other mechanism works, the ^ field 
settles down to this self-decaying configuration, whose energy density is however 
much more suppressed than that of the non-relativistic matter oc a~^. As will be 
discussed shortly, this does not mean that the particle creation permanently lasts, 
however. 

It is now appropriate to discuss how the resonant particle production ends and 
what the fate of the ^ field oscillation would be. As the ^ amphtude decreases both by 
cosmological damping and particle production, the 6 parameter given by 2/1^ /mj also 
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decreases, and when 6 becomes comparable to unity, relevant mode may enter into 
the small amplitude regime, or even into a stability band region. We suppose that the 
strong resonant decay effectively ends at this time and the ^ matter behaving like a 
non-relativistic matter persists until the Born decay lifetime becomes comparable to 
the Hubble time. Note in this context that within the first resonance band region the 
small amplitude decay rate formula, when translated per one particle by dividing 
the number density ra^ = | m^^^ , agrees with the Born formula yU^/ (327rm^). This 
means that the small amplitude region is well approximated by the Born formula, 
although the intermediate amplitude region of 6' ~ 1 is rather difficult to deal with. 

Let then the end time of resonant decay be denoted by r = . After this time and 
before the Born decay starts, the ^ matter density behaves like p^/pr oc a{T)/a{Te) . 
The Born decay starts around t = Td = 167rm|//i^ , from which one may estimate 
the fraction of ^ matter density at the epoch of Born decay, 

f = 16v/?/'/'^dV^ ^ = (|)"' f'^ . (20) 

Pr P^O <J Co "^pl 

with / = dTe of order unity numerically. This gives an upper bound on the pa- 
rameter combination, ■m^/{p'^^Q^'^m^(^) = I {gQrn^m}J^) , in order not to have a 
C dominant epoch. In the rest of discussion we assume that this bound is obeyed. 
Namely, we are considering a situation in which the ^ dissipation initially takes place 
via the rapid resonant particle production and at a later epoch a small ^ residual 
completely disappears via the Born decay. Details of this mechanism, in particular 
how individual Fourier modes undergo the Born phase of decay are left to further 
work. In actual situations there may well exist further strong dissipation of the ^ field 
due to the intermediate amplitude oscillation around 6 = 1. Here we simply assume 
that there is no sizable late time decay and ignore a possible spectrum component 
due to the Born decay. 

Details of the particle energy spectrum are obtained by solving the Boltzmann 
equation. A Boltzmann-like equation in our model of the expanding universe may 
be set up for the distribution function niu , t) of produced particles, once the source 
term is given by the reduced density matrix; 

d (X d 

u—n{uj,t)--uj^ — n{uj,t)=ujS[^], (21) 

with S[^] = ^tr (N^pd) = -j^Hi 2/p2i,2Z- The distribution function n{ijj ,t) here, 
when integrated over the entire phase space, is not normalized to unity as usual. It 
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rather gives the number density after the phase space integration. Solution to the 
Boltzmann-hke equation is then 

n(..0 = X'*'s(-||.*'). (22) 

from which one can compute the spectrum. 

Result of numerical computation is presented in a few figures. For lack of space 
we only show typical parameter cases of B/d — 10^ , 1 , 10~^ , and c^fg^jd — 10 , 
which we call large, intermediate, and small damping cases, all taking a large initial 
amplitude, = 200. As the initial condition to the differential equation we take for 
simplicity the ^ dominance so that the initial radiation energy density Pr(^o) = with 
the Hubble rate i?(io) = ■ The time scale in these figures is l/{dm^) — . 
We take this time scale larger than 1/m^ so that the dissipation time scale is much 
larger than the field oscillation period, for consistency of our approximation. In Fig.l 
is shown time variation of the ^ and the net radiation energy density, , Pr-^r^ 



a{t) 



along with the scale factor for the intermediate damping. Note that the asymp- 
totic behavior is observed already at moderately early times of a few times l/{dm^). 
In Fig. 2 is shown time variation of the energy spectrum in unit of m|/(27r^) plotted 
against the energy in unit of \/ h^q, in which is marked by the solid line the spectrum 
at the exit time of resonant particle production, taken for a guide to occur when 
9 = 1^ giving dr ~ 1.2 for this parameter set. The peak position of the spectrum 
shifts towards the low energy side as time proceeds, but the maximal peaking occurs 
at a time slightly earlier than the exit time. In Fig.3 we show the energy spectrum 
at the exit time in the three different cases of damping. The total number of created 
particles increases as the damping factor B/d decreases. 

Although the number distribution at the exit time is not the thermal one, there is 
a broad range of energies of which particles are produced, especially for larger initial 
amplitudes and larger damping. This tends to accelerate the approach to the ther- 
mal distribution in subsequent energy exchange process by interaction among created 
particles. The broadness of energy distribution is more pronounced for larger damp- 
ing or larger initial amplitudes, as illustrated in Fig.3. What is important to physical 
processes that occur at later times is the average energy and energy dispersion. We 
found numerically that in the dimensionless unit the energy dispersion approaches 
1, (o;^) — (o;)^ / {uj) — > 1 in the small B limit. This is a good measure of the 
broadness of the energy distribution. 
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A quantitative measure of dissipation is given by the entropy. Definition of the 
entropy in a classical system is somewhat ambiguous. In our approach based on 
the quantum behavior of radiation system it is natural to consider the quantum 
statistical entropy determined by the reduced density matrix p^- One may thus 
define the entropy of created particles by V" = — Tr {pu hi po) ) again the big 
trace including the phase space integral. In general it is not easy to derive a tractable 
expression of the entropy from this general formula, but use of the asymptotic form 
of the diagonal density matrix element for large /, P2i,2« ~ ^/^:^e~^', yields a 
simple, but suggestive form of s^^ ~ hi{Nj) for a single mode, ignoring a constant 
term of order 0.1. In terms of the previous variables, mode-summed entropy density 
at time t is then 

^ 1^""?^""'^^^ £ dT'Y{T')a\T'). (23) 

The entropy computed according to this formula is plotted in Fig. 4. It is seen that 
after initial transient behavior ends, the entropy decreases like ar^{t) with the scale 
factor, as expected. The asymptotic values of the net entropy s{t)a^{t) / a^{tQ) are 
4.0 , 43 , 110 in unit of m| for the large, intermediate, and small damping. 

For cosmological application it is important to be able to determine the reheat 
temperature after inflation. Even for a broad energy distribution of produced parti- 
cles, the energy redistribution towards the thermal one may be established only after 
a long elapse of time due to a small coupling. Nevertheless it might be reasonable, 
as a crude approximation in the case of broad particle distribution, to estimate the 
reheat temperature Tr by equating the radiation energy density at the exit time to 
1^ , the usual thermal energy density at the temperature Tr. This estimate gives 

T« = (^)^/^^^/^(re)yC^= ^(^)^/^o^/^5-V^ZV^(re)m,, (24) 



which in cases of most interest yields a fraction times \Jiom^. A similar estimate 
is obtained by assuming instantaneous thermalization and equating the Hubble rate 
H = \/ ^^o^ ~^ to the inverse of the exit time l/(2te), again giving Tr = a fraction 



times y^otn^. The reheat temperature thus derived is larger by a big margin than 
that given by the Born decay rate. The enhancement factor is roughly of order 
O[10] ^ '\J^i ■ ■^'-'^ instance, in a chaotic type of inflation with p = 10~^ , = 
10~^ rripi , ^0 — 0[mpi] , this factor may become as large as 10^, and Tr ~ ^J^om^ — 
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0[ 10^^ ] GeV . But in this particular example the energy spectrum is peaked around 
a low value of u/^/J^^ . Indeed, the concept of the reheating temperature is dubious 



The example of the chaotic inflation exceeds this limit. What happens for very small 
B is that the total energy density is not a good measure of the reheat temperature, 
because an enormous amount of particles is produced to give a large energy density. 
The distribution itself is broad, as evidenced by the unit energy dispersion and not 
too small an average energy. 

With a larger reheat temperature or more precisely with copious production of 
high energy particles of GUT scale energies, baryo-genesis may proceed efficiently in 
the original GUT scenario of asymmetry generation. A new serious problem would 
be a copious production of gravitino in supergravity models. The usual constraint 
on the reheat temperature favored by low energy supersymmetry is not met here 
due to the large number density of created particles. Resolution of this problem is 
left to future. 

How about the Polonyi problem [Q? We must admit that our present work does 
not deal with the intermediate case of ~ 1. Without considering a possible strong 
dissipation in the intermediate amplitude region the Polonyi field is damped by the 
factor -jl- 1^ j-^^g ^Yie resonant particle production. This is presumably 



not sufficient to solve the Polonyi problem, but we caution that the intermediate 
amplitude region must be included for a definite conclusion. 

In summary, we showed in this paper that cosmic field oscillation such as infia- 
ton or the Polonyi field is dissipated via rapid resonant particle production if the 
initial amplitude or the coupling to matter is large enough. The energy spectrum 
of produced particles along with detailed evolution of the field energy density may 
be followed by solving a closed set of differential equations in the large amplitude 
limit. Much larger production of high energy particles becomes possible than what 
the usual Born decay would imply. This may have great impact on baryo-genesis 
and gravitino production. 



for cases in which Tr is much larger than the kinematical limit of v^/i^o 
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Figure caption 



Fig.l 

Time evolution of ^ energy density (dash-dotted line) and radiation energy den- 
sity times scale factor to the 4th power (solid line) together with the Hubble rate 
a/ a (broken) for the intermediate set of parameters of c^/gQ/d = 10 ,B/d = 1 . The 
time scale in this and Fig. 4 is . The vertical scale is drawn in arbitrary 

units, with the ^ energy density given by the logarithmic scale on the right. 

Fig.2 

Time evolution of the energy spectrum in unit of m|/(27r^) plotted against the 
energy in unit of \//x|o for the intermediate damping case, as in Fig.l. The solid 
line is the spectrum at the exit time around rd — 1.2 , and the peak position shifts 
towards the low energy side as time proceeds according to rd — 0.5 , 0.75 , 1 , 1.5 , 2 , 
as indicated by dash-dotted and dotted lines. 

Fig.3 

Energy spectrum in unit of m|/(27r^) at various damping factors of B/d = 
10^ (dash - dotted) , B/d = l(solid) reduced by 10^ and B/d = lO-^(broken) 
reduced by 10^, plotted against the energy in unit of V/^Co ■ 

Fig.4 

Time evolution of the net entropy, entropy density times {t) / {tg) , for the 
same damping factor as in Fig.3. 
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